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Abstract 

We study the quasi-normal modes (QNM) of electromagnetic and 
gravitational perturbations of a Schwarzschild black hole in an asymp- 
totically Anti-de Sitter (AdS) spacetime. Some of the electromagnetic 
modes do not oscillate, they only decay, since they have pure imag- 
inary frequencies. The gravitational modes show peculiar features: 
the odd and even gravitational perturbations no longer have the same 
characteristic quasinormal frequencies. There is a special mode for odd 
perturbations whose behavior differs completely from the usual one in 
scalar and electromagnetic perturbation in an AdS spacetime, but has 
a similar behavior to the Schwarzschild black hole in an asymptotically 
flat spacetime: the imaginary part of the frequency goes as where 
r + is the horizon radius. We also investigate the small black hole limit 
showing that the imaginary part of the frequency goes as r\. These 
results are important to the AdS/CFT conjecture since according to 
it the QNMs describe the approach to equilibrium in the conformal 
field theory. 
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1 Introduction 



QNMs of black holes play an important role in the study of the dynamics 
outside black holes. They appear, for instance, when one deals with the 
evolution of some field in the black hole spacetime, or in black hole-black 
hole collision processes. Numerical simulations ranging from the formation 
of a black hole in a gravitational collapse QXj to the collision of two black holes 
[0] provide clear evidence that no matter how one perturbs a black hole, its 
response will be dominated by the QNMs. QNMs allow us not only to test 
the stability of the event horizon against small perturbations, but also to 
probe the black hole mass, electric charge and angular momentum, through 
their characteristic waveform. 

A great deal of effort has been spent to calculate the QNMs and their as- 
sociated frequencies. New powerful methods, both analytical and numerical 
have been developed. The main interest in these studies is in the application 
to the analysis of the data from the gravitational waves to be detected by the 
forthcoming gravitational wave detectors. We refer the reader to |3|, ffl for 
reviews. In a different context, York H tried to explain the thermal quan- 
tum radiance of a Schwarzschild black hole in terms of quantum zero-point 
fluctuations of zero mean in the QNMs. 

All these previous works deal with asymptotically flat spacetimes, but 
the recent AdS/CFT correspondence conjecture makes the investigation of 
QNMs in anti-de Sitter spacetimes more appealing. According to it, the black 
hole corresponds to a thermal state in the conformal field theory, and the 
decay of the test field in the black hole spacetime, corresponds to the decay 
of the perturbed state in the CFT. The dynamical timescale for the return 
to thermal equilibrium is very hard to compute directly, but can be done 
relatively easily using the AdS / CFT correspondence. Horowitz and Hubeny 
[0] (see also ||) began the study of QNMs in AdS, by thoroughly investigating 
scalar perturbations in 4, 5 and 7 spacetime dimensions. Subsequently, Wang 
and et al |9|, [Hj analyzed scalar QNMs in a Reissner-Nordstom AdS geometry. 
Recently, Cardoso and Lemos [jllj found an exact solution for the QNMs of 
scalar, electromagnetic and Weyl perturbations of a BTZ black hole. Another 
conjecture is related to the speculation 0, [12], [14]] that there might be a 
connection between the critical exponent of Choptuik [|13j and the imaginary 
part of the frequency, for small black holes. This is still an open question. 
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In this paper we shall go beyond the scalar perturbations |7|, ||, [HJ], and 
consider electromagnetic and gravitational perturbations of a Schwarzschild 
black hole in an asymptotically AdS spacetime. Electromagnetic pertur- 
bations are of interest due to the AdS/CFT conjecture since they can be 
seen as perturbations for some generic supergravity gauge field. In addition, 
the Maxwell field is an important field with different features from scalar 
or gravitational fields, which makes it worth studying. On the other hand, 
gravitational perturbations have the additional interest of arising from any 
other type of perturbation, be it scalar, electromagnetic, Weyl, etc., which in 
turn disturb the background geometry. Therefore, questions like the stability 
of spacetime for scalar or other perturbations, have a direct dependence on 
the stability to gravitational perturbations. 

We will find that in the case of electromagnetic perturbations of large 
black holes, the characteristic QNM frequencies have only an imaginary part, 
and scale with the horizon radius. As for gravitational perturbations, there 
are two novel features. First, contrary to the asymptotically flat spacetime 
case, odd and even perturbations no longer have the same spectra, although 
in certain limits one can still prove that the frequencies are almost the same. 
The second intriguing result is that, for odd perturbations, there is a mode 
with a totally different behavior from that found in the scalar and electro- 
magnetic case: in this mode the frequency scales with just as in asymp- 
totically flat Schwarzschild spacetime. We also investigate the small black 
hole limit (a problem recently addressed by Zhu et al |15|), and find that the 
QNM frequencies go as r\. 



2 Electromagnetic and Gravitational pertur- 
bations in a Schwarzschild AdS background 

2.1 Maxwell perturbations 

We consider the evolution of a Maxwell field in a Schwarzschild-anti-de Sitter 
spacetime with metric given by 

dr^ 

ds 2 = f(r)dt 2 - 9- - r 2 (d8 2 + sin 2 Odcf 2 ) , (1) 
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where, f(r) = (-^ + 1 — ^r), -R is the AdS radius and M the black hole 
mass. The evolution is governed by Maxwell's equations: 

F^. u = ,F ia , = A Vtll -Aw, (2) 

where a comma stands for ordinary derivative and a semi-colon for covariant 
derivative. As the background is spherically symmetric, we can expand A^ 
in 4-dimensional vector sphericall harmonics (see ||16|| ): 



E 

l.m 



V 



o 
o 

2 lm (t,r) 



l 9(pYi m 
a lm (t,r) sin 9d e Y l 



in i 
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f m (t,r)Y lm 
h lm (t,r)Y lm 
k lm (t,r)d e Y lm 
k lm (t,r)d^Y lm 



(3) 



/ 



where the first term in the right-hand side has parity (— and the second 
term has parity (—1)', vn is the azimuthal number and I the angular quantum 
number. If we put this expansion into Maxwell's equations (0) we get a second 
order differential equation for the perturbation: 



d 2 m 



drt 



+ 



- V{r) tf(r) =0 



(4) 



where the wavefunction \l/(r) is a linear combination of the functions / 



Im 



h , k and a as appearing in (Q). \I/ has a different functional de- 
pendence according to the parity: for odd parity, i.e, (— \P is ex- 
plicitly given by \I/ = a lm whereas for even parity (—1)' it is given by 
= (—iujh lm — ^-J , see |L6| for further details. It is assumed that 

the time dependence is $?(t,r) = e~ tuJt ^/(r). The potential V appearing in 
equation (|J) is given by 

~i(i + iy 



V(r] 

and the tortoise coordinate r* 



f(r) 



(5) 



is defined as 
Or 
dr* 

> and if we do this, the wave equation again 



f(r) 



(6) 



We can of course rescale r, r 

takes the form with rescaled constants i.e., r 
r + is the horizon radius. So, we can take R = 1 and measure everything in 
terms of R. 



oj — > W-R, where 
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2.2 Gravitational perturbations 



When dealing with first order gravitational perturbations one supposes that, 
at least in some restricted region of spacetime, the metric functions can be 
written as 

9ab{x u ) = g a ° b \x u ) + h ab (x u ) , (7) 

where the metric g a ° b \x u ) is the background metric, given by some known 
solution of Einstein's equations, and h a b{x v ) is a small perturbation. Our 
background metric is a Schwarszchild-anti-de Sitter metric ([!]) and the met- 
ric g a b{x u ) will follow Einstein's equations in vacuum with a cosmological 
constant: 

G ab ~ Agab = . (8) 

Upon substituting (0) in (P) we will obtain some differential equations for the 
perturbations. We use the same perturbations as originally given by Regge 
and Wheeler [18], retaining their notation. After a decomposition in tensorial 





19 and Mathews 



spherical harmonics (see Zerilli 
distinct classes - odd and even - with parities 



these fall into two 
{— and (—1)' respectively, 
where I is the angular momentum of the particular mode. While working 
in general relativity one has some gauge freedom in choosing the elements 
h a b{x u ) and one should take advantage of that freedom in order to simplify 
the rather lengthy calculations involved in computing ([8]). We shall therefore 
work with the classical Regge- Wheeler gauge in which the canonical form for 
the perturbations is (see also pi] ): 



odd parity: 



h 



[IV 







h (r) 







hi(r) 



h (r) 
h x {r) 


hoir) 



-iuit 



sin0^Wcos0); (9) 



even parity: 



H (r)f{r) 






#i(r) 
H 2 (r)/f(r) 








r 2 K(r) 









sin 



e' luJt Pi(cos6).(W) 
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Here Pi(cos6) is the Legendre polynomial with angular momentum /. If we 
put this decomposition into Einstein's equations we get ten coupled second 
order differential equations that fully describe the perturbations: three equa- 
tions for odd perturbations and seven for even perturbations. It is however 
possible to circumvent the task of solving these coupled equations. Regge 
and Wheeler [18] and Zerilli [23| showed how to combine these ten equations 
into two second order differential equations, one for each parity. So following 
Regge and Wheeler |L8| (see also |22| for more details) we define, for odd 
parity the wave function Q(r) given by perturbations, 



Q{r) 



fir) 



ht{r). 



(11) 



After some work, Einstein's equations yield 

d 2 Q 



dri 



+ 



w -V odd (r) Q = 



where 



V, 



odd 



f(r) 



1(1 + 1) 



Qm 



(12) 



(13) 



Likewise, following Zerilli |23| one can define for even modes the wavefunction 
T(r) implicitly in terms of Hq, Hi and K, through the equations 



K 



Qm 2 + c (1 + c) r 2 + m (3cr — 3 



2±\ 
b? ) 



r 2 (3m + cr) 
-3m 2 — 3cmr + cr 2 



3m w) 



r (3m + cr) f(r) 



T — iuj 



ar„ 

r dT 
f(r) dr* 



(14) 
(15) 



where c = | [1(1 + 1) — 2]. Then Einstein's equations for even parity pertur- 



bations can be written as 

o 2 t 



+ 



- V eveQ (r) 



T = 



(16) 



with 



2f(r) 9m 3 + 3c 2 mr 2 + c 2 (1 + c) r 3 + 3m 2 (3cr + 3-^r 



(3m + cr) 



(17) 
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Now, by defining 



«r 2771 ~ 3 
W = — + 



3r 



2c 3c 2 + 2c 2 + 27^ 

1 — + 3 

3c (3m + cr) 



(18) 



where j = -± + g + f^), we obtain 



odd 



dW 



dW 

dr* 



+ {3. 



(19) 



where (3 



9m 2 



It is interesting to note that the two potentials, odd 



and even, can be written in such a simple form, a fact which seems to have 
been discovered by Chandrasekhar ]2"5|] . Potentials related in this manner 
are sometimes called super-partner potentials ||24|| ). We note that similar 
equations were obtained by Mellor and Moss [^(| for Schwarzschild-de Sitter 
spacetime, using a different approach. 



3 Quasinormal modes and some of its prop- 
erties 

3.1 Analytical properties 

To solve equation (f|) for Maxwell fields and equations ([T2|-[n]) for gravita- 
tional fields, one must specify boundary conditions. Consider first the case 
of a Schwarzschild black hole in an asymptotically flat spacetime (see, e.g., 
[0]). Since in this case the potential vanishes at both infinity and horizon, 
the two solutions near these points are plane waves of the type \I/ ~ e ±iwr * , 
where the r* coordinate in this case ranges from — oo to oo. Quasinormal 
modes are defined by the condition that at the horizon there are only ingoing 
waves, i.e., \I/hor ~ e ~ lujr * . Furthermore, one does not want to have fields 
coming in from infinity (where the potential in this case vanishes). So, there 
is only a purely outgoing wave at infinity, i.e., ~ e wr * . Only a discrete 
set of complex frequencies u meet these requirements. 

Consider now a Schwarzschild black hole in an asymptotically AdS space- 
time. The boundary condition at the horizon is the same, we want that near 
the horizon \l/hor ~ e~ tur * . However, r* has a finite range, so the second 
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boundary condition needs to be modified. There have been several papers 
discussing which boundary conditions one should impose at infinity in AdS 
spacetimes ([^7|, We shall require energy conservation and thus 

adopt the reflective boundary conditions at infinity |27j. This means that 
the wavefunction is zero at infinity. For a different boundary condition see 



We now show that the imaginary part of the frequency lj is negative, 
for waves satisfying these boundary conditions, provided the potential V is 
positive. The proof proceeds as for the scalar field perturbation case [[?]], 
although there are some steps we think are useful to display explicitly here. 
Writing for a generic wavefunction as 

<f) = e iuir *Z, (20) 

where, Z can be Q or T, we find 

/W + [/ '-2 M f-^ = O, (21) 
where / = (r 2 + 1 — — ). In the proof, we are going to need the asymp- 



totic behavior of the solutions of equation fl2l|) . For r — > r + we have 
/ ~ (3r + + ^)(r — r + ) and X ~ C, where C is a constant which takes 
different values depending on the case, electromagnetic, odd or even gravita- 
tional perturbations. So equation (pip becomes, in this limit, 

Ay^ + [A-2i^-C^0, (22) 

where y = r — r + , and A = 3r + + — . This equation has an exact solution in 
terms of the modified Bessel functions I v (z) 131)1 , 



<V»I_2 (2(^2/)5) + Cay 4 */^ (2(^1/)*) • (23) 



We want the asymptotic behavior of these functions when y — > which is 
given by Ij, 
behaves as 



given by I v (z) — > r( y +1 ) , -2 — > 0. So, near the horizon the wavefunction <p 



(0) A A 

r (1 - ?f) r (1 + 
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We can see that if one wants to rule out outgoing modes at the horizon, we 
must have Ci = 0, so that in equation fl20|) does not depend on y. Let's 
now investigate the asymptotic behavior at infinity. For r — > 00 we have 
■7 — ► K^ 1 ) . Therefore near infinity equation ([21] ) becomes 



,<9 2 



+ [2r - 2iu;]-^- - 



Z(Z + 1) 



dr 2 1 1 dr 



0. 



Putting x = - we have 



5 2 



2 +2^^ -/(/ + 1)0 = 



with solution 



4>oo(x) = Ae 



dx 2 dx 
given by 

-*(m)) 1/2 



-lOJ — l\ CO 



-1(1+1)) 



1/2 



Now, 0oo(a; = 0) = 0, therefore A = —B, and thus, 



0oo(x) = Ae" 



sin 



;r 



-"(/«) 



Integrating by parts yields 



POO 

Jr+ 



dr 



dr 



dr 



~f\ 



— 2iuxf) 



dr 



dr 



— liujc 



7* 



d0 _ F 
dr / 



0. 







(25) 



(26) 



(27) 



(28) 



We can now proceed in the proof. Multiplying equation (|2l| ) by <fi (the 
complex conjugate of 0), and integrating from r + to 00 we obtain 



(29) 



(30) 



Now, one can show that ]$f-±\ r+ = 0, in order to satisfy the boundary 
conditions. Indeed, at r + , <f>(r + ) = constant and f(r + ) = 0. Now, at infinity, 
even though 0(oo) = 0, we have also /(oo) = 00, so we have to show that 
[4>fd]oa — 0. From equation (|28|) we can check that this is indeed true. 
Thus, equation ([$(]) gives 



dr 



dr 



dr } 



0. 



(31) 
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Taking the imaginary part of ([JT]) we have 



dr 



r+ 



UJ(f>— + Ul(j)—— 

dr dr 



(32) 



wich, after an integration by parts reduces to 



POO 




/ dr 


<j) — 


lr+ 


dr 



(33) 



Finally, inserting this back into fl3"T| ) yields 



dr 



Imw 



(34) 



From this relation, one can infer that, if V is positive definite then Imw < 
necessarily. So, since electromagnetic and even gravitational perturbations 
have V > one always has Imcu < 0. As for odd gravitational perturbations 
there are instances where V < 0, making this theorem unreliable for these 



cases. However, for r + < 



3 



, i.e., small enough masses, V > (see 
equation (0)), and the theorem applies. 

Another important point concerns the late time behavior of these fields, 



and the existence or not of power-law tails. As shown by Ching et al ||32|| , for 
potentials that vanish exponentially near the horizon, there are no power-law 
tails, so there will be no such tails in our case. 



3.2 Numerical Calculation of the QNM frequencies 

To find the frequencies u that satisfy the boundary conditions we first note 
that equation ([H]) has only regular singularities in the range of interest. It 
has therefore, by Fuchs theorem, a polynomial solution | ]33|| . To deal with 
the point at infinity, we first change the independent variable r to x = -. 
Now we can use Frobenius method by looking for an indicial equation (for 
further details see 0), and force it to obey the boundary condition at the 
horizon (x = ^ = h). We get the following solution to equation (pip, 

oo 

= Yj a «H ( x ~ h T ' ( 35 ) 

n=0 
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where a n ^) is a function of the frequency. If we put fl3~5|) into (^Tj) and use 
the boundary condition = at infinity (x = 0) we obtain 

oo 

IXhHO^ - ( 36 ) 

n=0 

Our problem is reduced to that of finding a numerical solution of the poly- 
nomial equation (|36|). The numerical roots for u of equation fl36|) can be 
evaluated resorting to numerical computation. Obviously, one cannot deter- 



mine the full sum in expression fl36|), so we have to determine a partial sum 
from to N, say and find the roots uo of the resulting polynomial expression. 
We then move onto the next term N + 1 and determine the roots. If the 
method is reliable, the roots should converge. We stop our search when we 
have a 3 decimal digit precision. 

3.2.1 Electromagnetic modes 

As long as the modes are decaying, it does not matter whether they're os- 
cillating or not. However, as we will see there are frequencies in the elec- 
tromagnetic case with a vanishing real part, which makes it possible to use 
an approximation, due to Liu [fL7|1 , to the highly damped modes. Although 
the method was originally developed for the asymptotically flat space, it is 
quite straightforward to apply it to our case. There is therefore a way to 
test our results. Unfortunately, this method relies heavily on having a fre- 
quency with a large pure imaginary part, so as we shall see it will only work 
for electromagnetic perturbations. We have computed the lowest frequencies 
for some values of the horizon radius r+, and I. The frequency is written 
as uj = uj r + iu>i, where uj r is the real part of the frequency and u>i is its 
imaginary part. In tables 1 and 2 we list the numerical values of the lowest 
quasinormal frequencies of electromagnetic perturbations for I = 1 and I = 2 
and for selected values of r + . For frequencies with no real part, we list the 
values obtained in Liu's aproximation. 
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Numerical 


Liu's approximation 


r+ 










0.8 


1.287 


2.175 






1 


1.699 


2.163 






5 


8.795 


~ 


7.6 


~ 


10 


15.506 


~ 


15.05 


~ 


50 


75.096 


~ 


75.01 


~ 


100 


150.048 


~ 


150.005 


~ 



Table 1. Lowest QNM of electromagnetic perturbations for I = 1. The — in 
Liu's approximation columns means that the method is not applicable. 





Numerical 


Liu's approximation 




-Ui 


UJ r 






0.8 


1.176 


2.501 






1 


1.579 


2.496 






5 


10.309 


0.822 


7.6 


~ 


10 


15.755 


~ 


15.05 


~ 


50 


75.139 


~ 


75.01 


~ 


100 


150.069 


~ 


150.005 


~ 



Table 2. Lowest QNM of electromagnetic perturbations for I = 2. 

As one can see, the imaginary part of the frequency, which determines how 
damped the mode is, and which according to the AdS/CFT conjecture is a 
measure of the characteristic time r = — of approach to thermal equilibrium, 
scales linearly (for large black holes) with the horizon radius supporting the 
arguments given in 0. Moreover, the frequencies do not seem to depend 
on the angular quantum number I, and are in excellent agreement with the 
analytical approximation for strongly damped modes. 
For a better visualization we also plot uj{ x r + in Figure 1. 
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Figure 1. Lowest electromagnetic QNM for / = 1 as function of r+. The 
black lozenges represent some frequencies numerically calculated. The line 
connecting them is a linear fit. 

3.2.2 Gravitational modes 

The numerical calculation of the quasinormal frequencies for gravitational 
perturbations proceeds as outlined previously (the associated differential 
equation has only regular singularities, so it is possible to use an expansion 
such as 

(i) Odd modes: In tables 3 and 4 we show the two lowest QNM frequencies 
for I = 2 and I = 3 odd gravitational perturbations. An important point 
in odd QNMs is that there is a mode for which the frequencies are not only 
pure imaginary and very small, but also scale with This is similar to the 
behavior of Schwarzschild black holes in asymptotically flat spacetimes, as 
mentioned. However, all frequencies have a negative imaginary part, which 
indicates that the spacetime is stable for this kind of perturbations. 
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lowest QNM 


second lowest QNM 


r+ 


-UJi 


UJ r 


-UJi 


UJ r 


0.5 


6.4 





0.72 


3.037 


1 


~2(?) 





2.404 


3.033 


2 


0.728 


~ 


5.258 


4.447 


5 


0.2703 


~ 


13.294 


9.577 


10 


0.13378 


~ 


26.626 


18.662 


50 


0.02667 


~ 


133.19 


92.505 


100 


0.0132 


~ 


266.384 


184.959 



Table 3. Lowest QNM of gravitational odd perturbations for 1 = 2. 





lowest QNM 


second lowest QNM 




-UJi 


UJ r 


-UJi 


0J r 


1 


10 





1.639 


3.849 


2 


2.189 





5.080 


4.615 


5 


0.690 


~ 


13.247 


9.735 


10 


0.336 


~ 


26.603 


18.742 


50 


0.0669 


~ 


133.19 


92.521 


100 


0.0333 


~ 


266.382 


184.967 



Table 4. Lowest QNM of gravitational odd perturbations for 1 = 3. 

The value u = 2 in Table 3 marked with a "?" is a somewhat dubious 



result. In fact, from (|24j) it follows that if 1 — = —n, with n an integer, 



then there is nothing going down the hole, so perhaps it is not a QNM. It is 
also an "algebraically special value" in the sense of Chandrasekhar [ |34|j . 

(ii) Even modes: In table 5 we show the lowest QNM frequencies for I = 2 
and / = 3 even gravitational perturbations. We point out the remarkable re- 
semblance of the values in table 3 with those in for scalar perturbations, 
even though the potentials are quite different. We have performed calcula- 
tions for higher values of the angular quantum number /, and found that the 
QNM frequencies are indeed very similar throughout all values of 



14 





lowest QNM, I = 2 


lowest QNM, I = 3 


r+ 






-UJi 




1 


1.584 


3.018 


1.392 


3.909 


2 


3.974 


4.546 


3.299 


4.597 


5 


12.649 


9.83 


11.642 


10.217 


10 


26.301 


18.806 


25.788 


19.089 


50 


133.125 


92.535 


133.022 


92.596 


100 


266.351 


184.959 


266.300 


185.005 



Table 5. Lowest QNM of gravitational even perturbations. 

(iii) Discussion: We first note that there is clearly a distinction between 
odd and even perturbations: they no longer have the same spectra, contrary 
to the asymptotically flat space case (see [0), a problem we shall consider in 
more detail in the next subsection. We also remark that in electromagnetic 
and scalar perturbations the frequency scales with r + (for large black holes 
at least). Since the temperature scales also with r + in the large black hole 
regime, this means that the frequency scales with the temperature. Thus, in 
the dual CFT the approach to thermal equilibrium is faster for higher tem- 
peratures. This is a totally different behavior from that of asymptotically 
flat space, in which the frequency scales with However, for odd modes 

there is one that scales with — . This is a reflection of the different behav- 

r+ 

ior of the potential V^d for odd perturbations (that was why we couldn't 
prove stability for odd perturbations in the first place), and of the boundary 
conditions, as we shall show in section 3.3. The odd modes are therefore 
particularly long-lived. 

3.3 On the isospectrality breaking between odd and 
even perturbations 

As is well known [p5| , |35| in the case of a Schwarzschild black hole in an 
asymptotically flat space the two potentials V even and V^d & ve T ^ se t° the 
same quasinormal frequencies (in fact to the same absolute value of the reflec- 
tion and transmission coefficients). This remarkable property followed from 
a special relation (the equivalent for asymptotically flat spacetimes of our 
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equation (|19D) between the potentials and the behavior of W at the bound- 
aries. However, as one can see in tables 3, 4 and 5 there is a isospectrality 
breaking between odd and even perturbations in Schwarzschild anti-de Sitter 
spacetime. 

We shall now treat this problem. The breaking of the isospectrality is 
intimately related to the behavior of W at infinity. On taking advantage of 
the machinery developed by Chandrasekhar, we seek a relation between odd 
and even perturbations of the form 

Q = PiT + Ql ^-, (37) 

T=p 2 Q + q2^-, (38) 



yielding (see |25[ for details), q\ = , Pi = qW , p 2 = —p\ and q 2 = qi 



q. Thus, we obtain 



dT 

Q = qWT + q—, (39) 

T = -qWQ + q^- . (40) 

Suppose now that u is a QNM frequency of T, i.e., one for which 

T - A eYen e-^ r * , r - r + , (41) 

T -> , r -> oo . (42) 



Substituting this into equation (pLOQ we see that 



Q ^ A cvcn g [W(r+) - iuj\ , r ^ r + , (43) 

'dT_\ 

dr * ),. 



Q^A^r) ,^oo. (44) 



However, from equation (P8|), f^-) is in general not zero so that will 
in general fail to be a QNM frequency for Q. Should Q and T be smooth 
functions of u, one expects that if q is "almost zero" then u should "almost" 
be a QNM frequency for Q. Now, the condition that q is almost zero is 
that (3 — uj 2 be very large, and one expects this to be true either when uj 
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is very large or else when f3 is very large. And in fact, as one can see in 
tables 3, 4 and 5 for very large u the frequencies are indeed almost identical. 
On the other hand, for very small black holes (j3 very large) one expects 
the frequencies to be exactly the same, since both potentials have the same 
asymptotic behavior in this regime, as we shall see in the next section. One 
would be tempted to account for the remarkable resemblances between QNM 
frequencies of scalar and gravitational perturbations by a similar approach, 
but the proof is still eluding us. Should such an approach work, it could 
be of great importance not only to this specific problem, but also to the 
more general problem of finding the asymptotic distribution of eigenvalues, 
by studying a different potential with (asymptotically) the same eigenvalues, 
but more easy to handle. 



4 The limit m 







Although it is not possible to solve exactly for the QNM frequencies, it is 
possible to gain some analytical insight in the special case of very small black 
holes. There has been some discussion about this regime (see [15|, and also 
[0j HH an d references therein). Here we shall exploit the behavior of QNM 
frequencies in this regime a little further. For very small black holes one can 
easily see that both potentials (electromagnetic and gravitational) look like, 
in the r* coordinate, a barrier with unequal heights: 



10+1) 



Figure 2. The potential for small black holes. V m 



4Z(l+l) 
27ri 



and a ~ 3r_i 



It is trivial to obtain equations for the quasinormal frequencies in this 
limit. If \I/ is a general wavefunction then 



— iu , r* < , 



(45) 
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d^/dr* %k x Be 



ik\Ce 



•qj Be ikir * + Ce~ ikir * 

d^jdr* ikDe ikr * - ik x Ee- ikr * 



De ikr * + Ee- 



ikr* 



< r* < a , 
, r* > a , 



(46) 
(47) 



where fci = (cu 2 — V max ) 1 ^ 2 , and fc = [cu 2 — 1(1 + 1)] . Imposing the continuity 
of the logarithmic derivative and furthermore that \& = at infinity (r* = |), 
we get 



k\+u) 

k\-w Jlikia I J 



e 2ifc(f-a) + 1 
1 _ e 2i*(5-o) 



(48) 



In the limit a — > 0, ki — > oo (m — > 0) we have, supposing that u> stays small, 
the condition e 2lk z = 1 which means that 



= An 2 + 1(1 + 1) 



n 



1,2, 



(49) 



corresponding to a bound state. This gives for the lowest QNM frequencies 
(n — 1): co>o = 2.45 for I = 1 and ujq = 3.16 for / = 2. The above are to be 
compared with those in Tables 1-5. The agreement seems excellent, and we 
can now go a step further : If we linearize ( |4*5| ) around the solution (^), i.e, 
if we write u = u + i5 and substitute back in (|8|) we obtain, to third order 
in <5, the values listed in Table 6 . 





a = 3/h 


a = Q/h 


1 


5 


5 


1 


-2.1/h-i(lA2/h' 2 ) 


-1.92/h-i(0.05/h 2 ) 


2 


-0.859/Zi - i(0.04//i 2 ) 


-0.85//i-i(1.4xl(T 4 //i 2 ) 



Table 6. The linearized frequency 5 for selected values of the angular quan- 
tum number / and the potential width a. 

We have chosen a typical value of a ~ |, but we can see that, although 
the real part does not depend very much on a, the imaginary part is strongly 
sensitive to a. Nevertheless, one can be sure that whatever value of a, the 
imaginary part goes as and is always negative. 
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5 Conclusions 



We have computed the electromagnetic and gravitational QNM frequencies 
of Schwarzschild-AdS black holes in four dimensions. These modes dictate 
the late time behavior of a minimally coupled electromagnetic field and of 
small gravitational perturbations, respectively. The conclusions are: (i) The 
frequencies all have a negative imaginary part, which means that the black 
hole is stable against these perturbations, since these will decay exponentially 
with time; (ii) Maxwell perturbations are strongly damped, so according to 
the AdS/CFT conjecture, any electromagnetic perturbed thermal state will 
rapidly approach equilibrium; (iii) for odd gravitational perturbations in the 
large black hole regime, the imaginary part of the frequency (decaying mode) 
goes to zero scaling with — , just as in asymptotically flat space. In terms 
of the AdS / CFT correspondence, this implies that the greater the mass, the 
more time it takes to approach equilibrium, an unusual result; (iv) scalar M 
and gravitational even perturbations exhibit an amazing similarity for the 
characteristic time damping of the perturbations, but we have not been able 
to prove it analytically; (v) in the small black hole regime the imaginary part 
of the frequency (decaying mode) scales with r+. 
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